Nonperiodic echoes from mushroom billiard hats 
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Mushroom billiards have the remarkable property to show one or more clear cut integrable islands 
in one or several chaotic seas, without any fractal boundaries. The islands correspond to orbits 
confined to the hats of the mushrooms, which they share with the chaotic orbits. It is thus interesting 
to ask how long a chaotic orbit will remain in the hat before returning to the stem. This question 
is equivalent to the inquiry about delay times for scattering from the hat of the mushroom into 
an opening where the stem should be. For fixed angular momentum we find that no more than 
three different delay times are possible. This induces striking nonperiodic structures in the delay 
times that may be of importance for mesoscopic devices and should be accessible to microwave 
experiments. 

PACS numbers: 05.45.Gg, 05.45.Pq 



I. INTRODUCTION 

Two-dimensional billiards play a central role in the de- 
velopment of chaos theory ever since the early work of 
Sinai [l[ . In physics they have acquired increasing impor- 
tance as they are seen to emulate properties of systems as 
different as quantum dots [H or planetary rings [§] . Par- 
ticularly quantum or wave realizations of such objects 
have become popular, since fiat microwave cavities, so- 
called microwave billiards, are used by experimentalists 
at different laboratories 0, [E B S H B El > and these ex- 
periments mimic some properties of mesoscopic devices. 
From a mathematical point of view one of the advantages 
of billiards is that, in many instances, chaotic properties 
can be proven for cases of complete chaos 
and more recently even for mixed systems 




Microwave billiards serve as an analog system for the 
experimental study of the wave behavior of the cor- 
responding classical billiard, and permit a direct test 
of hypotheses proposed concerning connections between 
the quantum and the corresponding classical dynamic 
[l6l fl7l ] . In particular the spectral behavior and prop- 
erties of the wave functions of quan tum s y stems, whose 
classical analogue is chaotic fT§l. Il9l l20l l2ll I22I l23j , inte- 
grable [24l . [2{| or intermediate |26l. l27l [281] can be inves- 
tigated experimentally with such systems. In scattering 
systems similar matters have been discussed mainly for 
systems that are chaotic and hyperbolic [H|,[2!|, or mixed 
[30l . l3ll . l32l [H, HH ■ The role of parabolic manifolds has 
also received considerable attention 0, [H| . 

Mixed systems typically have the property, that in 
some border region integrable and chaotic areas inter- 
mingle as a fractal. Such a scenario is well illustrated by 
the twist map [361 ] . It may help to detect some character- 
istics of the dynamics in certain cases [34j . Bunimovich 
[Tij has proposed a family of mushroom billiards with 
mixed phase space which have the unusual property that 
chaotic areas and integrable ones are not separated by a 
fractal set of integrable islands extending into the chaotic 



sea. Such billiards are characterized by circular or elliptic 
hats connected to a stem or stems composed of straight 
walls; we shall concentrate on the case with only one hat 
and one stem. Typical examples are shown in Fig. [TJ The 
stem pertains entirely to the chaotic area, while the hat 
houses both chaotic and integrable trajectories. Several 
studies on the classical dynamics in mushroom billiards 
have already been presented [13, [38| . 

It is clear, that the relevant properties of such a bil- 
liard are determined by the hat of the mushroom, and in 
view of the successful theoretical [33| and experimental 
[34j analysis of phase space structures in the context of 
scattering echoes we wish in the present paper to ana- 
lyze the hat of the mushroom alone, viewed as a scatter- 
ing system with an opening where the stem used to be; 
yet to establish contact with other work, we shall also 
keep the possibility of a stem in mind. In the spirit of 
scattering echoes we ask how long a particle stays in the 
hat, if injected from the opening of the billiard or from 
the stem. Note that in the case of a semicircle mush- 
room hat with a symmetric stem, the limit of the area of 
stable bounded orbits is a caustic, which forms a circle 
segment connecting the points where the stem starts, as 
seen in the mushroom billiard in Fig. QJa). The usual 
Smale horseshoe construction [l5| fails, because there is 
no hyperbolic periodic orbit along this line. We have 
therefore numerically investigated the distribution of the 
classical delay times of such a billiard and found a sur- 
prising selectivity in the allowed delay times. 

As the angular momentum is a constant of motion in 
the scattering process as long as a particle stays within 
the hat we shall first consider fixed angular momenta and 
we find that generically only three delay times occur. On 
a subset of measure zero, one and two delay times are also 
possible; a larger set of different delay times is outright 
forbidden. This result is intimately related to old results 
on the circle map [3!| and allows to understand the above 
mentioned selectivity. 
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FIG. 1: Different shapes of mushroom hats. Circular mush- 
room billiards consist of a hat with a circular boundary and a 
stem, which has been chosen in triangular form for all mush- 
room billiards here. Stable orbits in a circular mushroom hat 
never enter the stem. They form a caustic in the hat, and the 
radius of the smallest possible caustic is one-half the width 
of the stem, as indicated in (a). Other shapes of mushroom 
billiards may display an elliptic hat or a shifted stem (b), 
very wide or narrow stems (c) or hats which are a circle sec- 
tor of an angle smaller than 180° (d), in this case 90° and 
180° x (V5-l)/2. 



II. DEFINITIONS AND PRELIMINARY 
RESULTS 

The circular hats of mushrooms are characterized by 
two quantities: The position and size of the hole and 
the angle between the straight walls that constitute the 
underside of the hat. In Fig. [1] we show various hats 
of mushroom billiards with different properties: circular 
mushroom hats [(a) and (c)], mushroom billiards with 
a shifted stem or an elliptic hat (b) and mushroom bil- 
liards whose hat is a circle sector of an angle different 




FIG. 2: The O billiard is a circle billiard with a straight line 
of length 2r along a diameter of the circle of length 2R, which 
defines the starting and ending points of the particle orbits 
taken into consideration. In a semicircle mushroom hat with 
the opening located at the position of the straight line these 
orbits correspond to those of particles entering and leaving it. 
In this sketch, a particle starts from the straight line at Ps, 
undergoes four reflections with the boundary at P\, P2, P3, 
and P4, before it finally ends on the straight line at Pe- 



from 180° (d). The stem is always chosen in triangular 
form to avoid the trivial parabolic manifold of so-called 
" bouncing ball" orbits between parallel walls of a rectan- 
gular stem. 

In Fig. [TJc) we show 180°-mushroom billiards with a 
wide and a narrow stem, respectively. Here a narrow 
stem is one narrow enough, that orbits which are trian- 
gular in the full circle remain within the hat irrespective 
of the orientation of the triangle. This implies, that all 
polygons, which do not form stars with intersecting line 
segments are confined to the hat. Moreover, the delay 
time, i.e., the time a particle entering the mushroom hat 
stays there, is approximately determined by the number 
of its bounces off the circle boundary. This is clearly 
not the case for a wide stem, where some polygons are 
no longer confined to the hat. In the following we will 
restrict ourselves to mushroom billiards with a narrow 
stem and count bounces of particles entering and leaving 
the mushroom hat, instead of determining delay times. 

We also show mushroom billiards where the two 
straight walls constituting the underside of the hat have 
angles smaller than 180° in Fig. QId). Note that it is 
essential that the two straight walls are parts of radii of 
the circle, i.e., form an angle of 90° with the tangent to 
the circle boundary at their intersection. Otherwise we 
lose the properties of the mushroom billiard. If the an- 
gle is larger than 90°, such that the extensions of the 
two straight walls meet in a point below the center of 
the circle, the system becomes ergodic and chaotic [40| . 
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For an angle smaller than 90° we get a more complex 
coexistence of integrable and chaotic dynamics. 

In the following we first will consider open ISC- 
mushroom billiards with a semicircular hat and a sym- 
metric stem. We are interested in the number of bounces 
a particle entering a mushroom hat through the opening 
experiences at the circle boundary before exiting again. 
Bounces with the straight part of the boundary are not 
counted. Whenever the particle hits it, we may use the 
principle of mirror images (4lj |. i.e., reflect the semicircle 
across its straight boundary, thereby obtaining a com- 
plete circle and continue the orbit into its lower half. Ac- 
cordingly, we may introduce a convenient alternative sys- 
tem, the so-called 9 billiard [37| shown in Fig. ® which 
consists of a circle billiard, whose radius R equals that 
of the mushroom hat, and a straight line of length 2r 
where the opening of the corresponding semicircle mush- 
room hat is. The only purpose of the straight line is, 
that it defines the starting and the endpoints of those 
particle orbits we are interested in. In the mushroom 
billiard these correspond to particles entering and leav- 
ing the semicircle hat. Accordingly, in order to compute 
delay times, we consider particle orbits starting on the 
straight line in the O billiard and count the number of 
bounces of the particle with the circle boundary, until 
it again reaches the straight line. We shall see that the 
© billiard can be readily generalized to describe arbitrary 
mushroom hats. While the © billiard is convenient for 
mathematical purposes physically it seems of little direct 
use. 

As already noted above and indicated in Fig.flja), for 
every particle orbit in the circle billiard the line segments 
connecting subsequent reflection points form a caustic of 
circular shape around the center of the billiard. The 
radius of the caustic equals the distance of these line seg- 
ments from the circle center. Therefore it gives the angu- 
lar momentum of the particle where we set the absolute 
value of the momentum to unity, as it is a conserved 
quantity. If the angular momentum of the particle is 
large enough, the corresponding orbit will never hit the 
straight line in the billiard. Such orbits belong to the 
integrable part of the phase space of the mushroom bil- 
liard. If they are periodic they may be polygons or stars. 
If they are aperiodic they will correspond to slowly rotat- 
ing stars or polygons, i.e., they would be periodic in some 
rotating frame of reference. The limiting angular mo- 
mentum corresponds to the radius of that caustic, which 
just touches one or both ends of the straight line. Any 
aperiodic orbit with angular momentum smaller than the 
limiting value will eventually reach the straight line after 
some number of bounces. For periodic orbits we can have 
the particular case, where a star or polygon will intersect 
the straight line in certain orientations and not in oth- 
ers. The orientations that do not lead to intersections 
will yield a parabolic manifold in the chaotic sea as ex- 
plained in Ref. [37| • This paper focuses on orbits that 
do cross the straight line. 

We now come to a really surprising result. Using a 



simple reflection program, which computes orbits of a 
particle in a billiard based on the law of specular re- 
flection, we counted the bounces off the circle boundary 
a particle starting from the straight line with randomly 
chosen initial conditions experiences until it reaches the 
straight line again. We found a strange selectivity in 
these numbers, which resembled a generalized Fibonacci 
sequence. Indeed for r = R/'i the first observed numbers 
of bounces are 1, 4, 5, 9, 14, 23, 37. The next number in 
this sequence, 51, does not continue the Fibonacci series, 
but is the sum of 37 and 14. Similar results have been 
found for other ratios r/R. We shall call the possible 
numbers of bounces. 

The purpose of this paper is to understand the ob- 
served selectivity. For this we will introduce a map, which 
generates the dynamics in the mushroom hat in an effi- 
cient way and then shall proceed in two steps utilizing 
this map. First we shall analyze the problem for a fixed 
angular momentum, considering that angular momentum 
is a conserved quantity both in the open mushroom hat 
and in the billiard. We shall obtain the surprising re- 
sult that typically only three magic numbers termed a 
magic triplet are possible for each angular momentum. 
To understand the entire delay time structure, we then 
analyze how these triplets evolve when changing the an- 
gular momentum. Based on the map, we developed a 
very efficient algorithm for the computation of the magic 
numbers. 



III. MAGIC NUMBERS FOR FIXED ANGULAR 
MOMENTUM 

Let the radius of the billiard be R and the straight 
line of length 2r be oriented symmetrically along the hor- 
izontal axis. A particle orbit is characterized by the par- 
ticle's angular momentum M, which is just the radius of 
its caustic. Orbits with angular momentum M > r never 
intersect the straight line, whereas orbits with M < r 
intersect it (up to a set of marginally unstable orbits of 
measure zero as stated above) at some time and therefore 
correspond to chaotic trajectories in the mushroom bil- 
liard. The angular momentum can be expressed in terms 
of the angle between the orbit of the reflected particle and 
the tangent to the circle boundary at the point of impact, 
which also is a constant for each particle orbit; we will 
call this angle tp, where < tp < n/2 (see Fig. [3]). Then, 
the radius of the caustic, i.e., the angular momentum M 
equals 

M = Rcostp. (1) 

A particle orbit is also characterized by its initial orien- 
tation. Generally, the orientation of the particle orbit 
after n reflections at the circle boundary may be given 
in terms of the orbit's angle with the positive horizon- 
tal axis. Equivalently, its orientation may be defined in 
terms of the total angle covered by the vector pointing 
from the circle center to the point of contact between the 
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caustic and the orbit during its rotation while the par- 
ticle is being reflected at the circle boundary. For our 
purposes it is more convenient to define it in terms of 
the latter. We denote this angle by n for the nth line 
segment of a particle orbit, which connects the nth point 
of impact with the (n + l)th. At each reflection the angle 
increases by 2ip (see Fig. Then, for a particle orbit 
starting from the straight line with an initial orientation 
#o measured with respect to the positive horizontal axis 
and an initial angular momentum M, which defines ip, 
()„ reads 



■ n2ip. 



(2) 



must evaluate this inequality for all allowed initial ori- 
entations 6*o [see Eq. ([3])] and resulting orientations 8 n . 
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Note, that the angle 9 n is not restricted to the interval 
[0, 271"), as it also contains the information on the number 
of times the particle travelled around the circle center. 
In conclusion, each particle orbit is completely character- 
ized by the initial values (if),0o), the radius of the caustic 
[Eq. {T)] and the orientation of each of its line segments 
[Eq. ©]• 

Due to the symmetry of the system it is sufficient to 
deal with orbits proceeding from the straight line into the 
right upper quarter of the billiard. For a given initial 
orientation #o the orbit starts on the horizontal line at a 
distance Mj cos6>o from the circle center (see Fig. [3]). As 
we are only interested in orbits starting from the straight 
line, which reaches from — r to r, we obtain as a necessary 
condition for the initial angle 9q 



-x< 



< 



X 



with the notation 



X = arccos(M/r), 



where with Eq. (Q]), 



< X < i>- 



(3) 
(4) 
(5) 



FIG. 3: A 9 billiard with the straight line chosen symmetric 
with respect to the center of the circle. Initial conditions for 
a particle orbit in the Q billiard are the angular momentum 
M, it i.e., the radius of its caustic (dotted line circle), and 
the orientation 8o of its first line segment. The orientation of 
each line segment of a particle orbit is defined in terms of the 
total angle covered by the vector pointing from the center of 
the circle to the point of contact of the line segment with the 
caustic during its rotation while the particle propagates. In 
this example the orientation 8o of the initial line segment is 
negative and determined by extending the latter back beyond 
its point of contact with the caustic as indicated in the figure 
by the dotted line. The angle tp between the reflected particle 
orbit and the tangent to the circle boundary at the point of 
reflection is a constant. At each reflection the orientation 
angle increases by 2ip. 



For an angular momentum M and length 2r of the 
straight line, ±% in Eq. @ correspond to the initial ori- 
entation of those particle orbits which start at the end- 
points of the straight line. For a particle starting with 
an initial orientation Oq the nth line segment of its orbit 
intersects the horizontal line at Mj cos#„, if measured 
with respect to the circle center. Hence, it crosses the 
straight line, if — r < Mj cos#„ < r, i.e., if [with Eq. 

0„mod27re [0; X ] U [tt - x; + X] U [2tt - x; 2tt] , (6) 

or equivalently if there is an integer number m such that 

TO7T — x < n < mir + x , m — 0, 1, • ■ • , (7) 

that is, if the distance between 9 n and an integer mul- 
tiple of 7T is smaller than x- I n order to determine all 
possible numbers of bounces that particles starting from 
the straight line perform until they reach it again, we 



With Eq. Q and Eq. |2]) the latter take values from in- 
tervals of length 2x around integer multiples of 2-0, 

n2ip-x < n < n2ip + x , n = 0, l,---. (8) 

While this inequality defines the range of possible values 
for the orientations 6 n , that in Eq. ([7]) gives the condition 
for the intersection of a particle orbit with the straight 
line after n reflections. Whenever one of the intervals 
defined in Eq. ([5]) has common values with one of those 
given in Eq. (0, a part of all possible orbits reaches the 
straight line. Thus, the only information we need for the 
computation of the magic numbers are the orientations of 
the line segments of a particle orbit. This procedure for 
the computation of the magic numbers is fast and much 
more efficient than the reflection program. 
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FIG. 4: Illustrative representation of particle orbits starting from the straight line in the Q billiard in terms of the orientations 
of their line segments. The orientations of the ensemble of all particle orbits starting from the straight line are restricted to the 
interval [— x>x] [ see Eq. © in the main text]. Their evolution is illustrated by shifting this interval by 2ip at each reflection of 
the particle with the billiard boundary, thereby obtaining the range of values of the orientations of successive line segments of 
the particle orbits [see Eq. ([8jl]. Thus, a string of "orbit intervals" is obtained. The straight line is represented by a string of 
"line intervals" of size 2\ around each multiple of n on the angle axis. When an overlap of a line and an orbit interval occurs, 
the overlapping part of the orbit interval corresponds to angles 9 n , which fulfill the inequality Eq. iffjl. that is, to initial angles 
do, for which the particle orbit ends after n reflections. Accordingly, the overlapping part is cut off the orbit interval in the 
subsequent iterations. 



For a given angular momentum M, which defines the 
angle ip, and length of the straight line 2r giving the 
angle x, an illustrative graphical representation of this 
inequality is obtained as follows: All possible initial ori- 
entations 6q of the particle orbits we are interested in 
take a value from an interval of length 2\ situated sym- 
metrically around [see Eq. (J3J], According to Eq. ([8]) 
the orientation 6 n of the nth line segment of a particle 
orbit then takes a value from an interval of length 2\ 
shifted by n2ip with respect to the initial interval and by 
2ip with respect to that of the preceding line segment, 
6 n -i. Hence, the orientations of the line segments of all 
those particle orbits, which start from the straight line, 
resume values from a string of angle intervals of length 
2\ situated at distances of 2ip along the angle axis as 
depicted in Fig. [U We call these intervals orbit inter- 
vals. Next, let us draw a string of intervals of length 2\ 
at distances ir along the angle axis below the string of 
orbit intervals; we will call these intervals line intervals. 
While the orbit intervals refer to the inequality Eq. © 
the latter refer to the inequality Eq. ([7]). The inequality 
Eq. ([7]) is fulfilled for a part of the interval of allowed 
initial angles, whenever one of the orbit and one of the 
line intervals have common values, i.e., partially overlap. 
As we are only interested in orbits which start from the 
straight line, the first orbit interval completely overlaps 
with the first line interval. 

In the example shown in Fig. [?J an orbit interval over- 
laps with a line interval after the first shift of the orbit in- 
terval. Hence, a part of the incoming particles reaches the 



straight line after one bounce with the circle boundary 
and the corresponding orbits end there. Accordingly, the 
orbit interval is shortened by this part. After four reflec- 
tions, another part of all possible particle orbits reaches 
the straight line, i.e., the orbit interval is further short- 
ened. The remaining part of the orbit interval overlaps 
with a line interval after five reflections. Thus, in this 
example the magic triplet consists of the numbers 1, 4, 
and 5. 

The Eqs. ((2|) and J6]) are related to old results on the 
so-called circle map [39]. On the basis of these results it 
can be shown that, for a fixed angular momentum, there 
are at most three magic numbers, and that if there are 
three magic numbers the largest one is the sum of the two 
smaller ones. We found an alternative proof of this result 
which in fact generalizes the translation of orbit and line 
intervals as depicted in Fig. |4j Using this approach we 
also obtained the following results: There is a subset of 
measure zero, where we have only two magic numbers. 
Further, in the case of the parabolic manifolds described 
in Ref. (37j only one magic number is finite, whereas the 
others are infinite and correspond to marginally unstable 
orbits which never leave the mushroom hat. Finally, in 
the case of angular momentum zero, all particle orbits 
intersect the straight line after a single bounce. We thus 
have exhausted all possibilities. We give in Table U the 
magic triplets for different angular momenta for r = R/3 
used in Fig. [IJa), and we readily recognize the numbers 
therein. 



IV. MAGIC NUMBERS AS A FUNCTION OF 
ANGULAR MOMENTUM 

As the appearance of only one or two magic numbers 
is restricted to a subset of measure zero in the set of 



possible angular momenta we must vary the angular mo- 
mentum in order to investigate them. In Fig. [5] we show 
the magic numbers as a function of angular momentum 
for the semicircle mushroom hat with radius R and half- 
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TABLE I: Magic triplets for different values of angular momentum M for a circular mushroom hat with radius R and half- 
opening r = R/3. 



M 


0.1000 


0.5000 


0.9000 


0.9500 


0.9800 


0.9950 


0.9980 


0.9990 


0.9992 


ni 


1 


1 


1 


4 


5 


9 


14 


14 


37 


n> 


4 


4 


4 


5 


9 


14 


23 


37 


51 


713 


5 


5 


5 


9 


14 


23 


37 


51 


88 



opening r — R/3. We observe discontinuities in the be- 
havior of the triplets of time delays (some are marked 
with arrows in the figure). Precisely there one member 
of a magic triplet vanishes, such that there are just two 
magic numbers, and a new magic triplet emerges. The 
new magic number equals the sum or difference of the 
other two. This apparently explains the preliminary re- 
sult that within the whole spectrum of magic numbers 
each except the two smallest is the sum of two smaller 
ones. 

Moreover, we observe a singularity in Fig. [5j It is 
located at an angular momentum, which is associated 
with a family of marginally unstable periodic orbits [37| . 



1000 




0.33300 



0.33330 



Angular momentum / R 



FIG. 5: Magic numbers plotted as a function of angular mo- 
mentum for a semicircle mushroom hat with radius R and a 
half-opening r — R/3 chosen symmetrically with respect to 
the circle center. The angular momentum is given in units of 
R as the absolute momentum is set to unity. One observes 
lines where a magic number remains constant over a certain 
range of angular momentum values, and sudden discontinu- 
ities (arrows), where one member of the triplet changes its 
value. Note further a singularity, where two of the three mem- 
bers of the triplet tend to infinity. Singularities of this kind 
are observed at angular momenta of unstable periodic orbits, 
which do not intersect the straight line in the O billiard in 
some orientations and do it in others. Between singularities 
of this kind, every magic number except the first and second, 
is the sum of two smaller ones, as proven in the text. Note 
that the chosen interval of angular momentum is small and 
corresponds to particles which reach the straight line close to 
one of its ends. 



These periodic orbits are marginally unstable, because 
the smallest change in angular momentum will cause 
them to rotate about the circle center and thus eventually 
to hit the straight line in the 9 billiard. The smaller the 
change in angular momentum, the slower this rotation 
will be, and this must lead to the delay times diverging 
to infinity as we approach the angular momentum of the 
periodic orbit. This is clearly seen in Fig. [5] Note that 
the slopes of the flanks formed by the discontinuities to 
the left and the right of the angular momentum asso- 
ciated with the singularity are different. Between zero 
angular momentum and the first singularity as well as 
in each interval between two singularities only disconti- 
nuities of the type discussed above may happen. There- 
fore, within an interval, each magic number except the 
two smallest ones at the minimum equals the sum of two 
smaller magic numbers. However, there is no such rela- 
tion between magic numbers from two different intervals. 

The fact that most magic numbers can be decomposed 
into two smaller ones, immediately explains why we ob- 




Number of bounces 

FIG. 6: Counts of the numbers of bounce particles starting 
from the straight line of half-length r — R/3 in the bil- 
liard experience with the circle boundary until they reach the 
straight line again; they were obtained using the linear in- 
equality Eq. (0 for the orientation of line segments between 
two reflections. Note the persisting selectivity. The organized 
and regular long time behavior (large number of bounces) is 
dominated by unstable periodic orbits, while the short time 
behavior (small number of bounces) cannot be explained by 
those. 
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for all triplets with first magic number one, i.e., if 
M < ' 



y/l + 1/r 2 

Using the sum rule, the third magic number is given by 

m = 112 + 1. 

For very small angular momenta, the angles ip(M) and 
x(M) can be expanded in a Taylor series, leading with 
Eq. © to the simple expressions n\ — 1, n 2 — [1/r] + 1, 
and 713 = [1 /r] + 2 for the very first magic triplet within 
a spectrum of magic numbers of an arbitrary mushroom 
hat. 



FIG. 7: Counts of magic numbers for mushroom hats with 
half-opening r = 0.3R and different shapes: 180°-mushroom 
hat (open triangles), 180° x (\/E — l)/2-mushroom hat (open 
circles), and 90°-mushroom hat (open squares) [see inset and 
Figs. 1(a) and 1(d)]. Note that the selectivity is observed in 
all examples, as the density of points remains almost constant 
in the logarithmic plot. For wider openings there is no corre- 
spondence between the number of reflections and the physical 
time any more. 



serve a strong selectivity in the spectra of magic num- 
bers. In Fig. [6] we show the counts of bounce numbers for 
r = R/3 in a double logarithmic plot up to 50000 bounces 
obtained using Eq. (J7J). We clearly see that scarcity dom- 
inates the picture. Note that the sequence of the magic 
numbers starts with the Fibonacci like behavior as dis- 
cussed above. Large magic numbers are dominated by 
the marginally unstable periodic orbits, which cause two 
organized sequences of magic numbers with a certain pe- 
riodicity for each singularity, one for angular momentum 
values below, the other for values above the singularity. 
Therefore we understand, why magic numbers may con- 
sist of subsets which show a periodicity as observed in 
Ref. [33]. For our purpose, Fig. [5] teaches us two facts. 
First, that we are really able to solve the problem also 
for very large bounce numbers. Second, that the small 
bounce numbers display a selectivity, that may be ac- 
cessible to experiments as the allowed numbers typically 
translate into fairly distinct physical times. 

It turned out that some magic numbers of a mushroom 
hat can actually be calculated analytically. For every 
mushroom with a central stem the first magic number 
equals one, ni = 1, as in the limiting case of vanishing 
angular momentum, every orbit leaves the mushroom hat 
after one reflection. Furthermore, the second magic num- 
ber of all triplets, whose first magic number is n\ = 1 is 
given by 



where [x] denotes the integer part of x. This result is 
obtained by evaluating the inequality Eq. IjTj). It is valid 



V. GENERALIZATIONS TO OTHER 
MUSHROOM BILLIARDS 

All results presented so far were given for a 180°- 
mushroom billiard implying a period of 7r for the line 
intervals. For hats with an angle SI smaller than n the 
entire argumentation holds as above except that the pe- 
riod 7T is replaced by f2. Note that this corresponds to a 
" generalized 8 billiard" , where the straight line rotates 
each time the straight wall constituting the "underside" 
of the hat is hit. If the mushroom hat has an asymmet- 
ric opening, the straight line must be shifted even for the 
180°-billiard, as the asymmetry is reversed each time the 
particle hits its underside. The chain of arguments be- 
comes a little more tedious, but still goes through. Thus 
the description of the dynamics inside the mushroom hat 
in terms of the translation of orbit and line intervals is 
indeed general, such that the results hold for a much 
wider range of mushroom billiards. In Fig. [7] we show 
bounce numbers in a double logarithmic plot for several 
angles and openings. The important point is, that we see 
no qualitative difference in the short time behavior: the 
selectivity remains in all cases. Furthermore, the magic 
numbers again combine to triplets (not shown in the fig- 
ure). 

We may finally go one step further and ask, what hap- 
pens with elliptic hats [421 ]. There we must remember, 
that such a hat displays two families of stable particle 
orbits, one corresponding to particles travelling between 
the two focal points, the other one to particles surround- 
ing them. We performed numerical calculations for three 
different choices of the positions of the straight lines for 
an elliptic billiard using the reflection program. In the 
first case, the straight line was positioned between the 
two focal points and in the second it extends across the 
focal points. The results are shown in Fig. [8l We see 
no selectivity in the first case (circles) and a picture very 
similar to those obtained for the semicircle hats in the 
second case (crosses). In the second case particles mov- 
ing on orbits, which start between the focal points all 
reach the straight line again after only one bounce with 
the ellipse boundary. The orbits surrounding the two 
focal points on the other hand are just deformed stars. 
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FIG. 8: Magic numbers of an elliptic billiard, correspond- 
ing to a mushroom billiard with an elliptic hat [see Fig. 1(b)], 
with the length of the semimajor axis chosen equal to R, that 
of the semiminor axis equal to y/3/2R while the focal points 
have a distance of 0.5R from the center. Circles refer to a 
straight line chosen symmetrically between the focal points 
(length 2r = 0.99/?), where no selectivity is observed. For a 
straight line of length 2r = 1.1 R selectivity is evident even 
for long times, because here only deformed starlike periodic 
stable orbits of the ellipse survive. Those particles, which 
intersect the straight line between the focal points, will inter- 
sect it again after just one bounce with the elliptic boundary. 
The diamonds refer to the case of two straight lines of length 
0.1R, whose outer endpoints are located in the focal points. 
They cut the caustic of the stable orbits of particles crossing 
the horizontal line between the two focal points from outside. 
Again a selectivity is observed. 



The line segments of such stars and therefore their ellip- 
tic caustic cross the straight line in a way similar to that 
of those in a circular 0-billiard. This brought us to the 
assumption, that a necessary condition for the selectivity 
is that the opening of the mushroom hat cuts a limiting 
caustic from outside. In order to confirm this picture 
we performed a third numerical simulation, where two 
straight lines were positioned in the elliptic O billiard, 
which extend from the focal points to the interior (see 
inset in Fig. [5]). In this case the hyperbolic caustic of 
those stable orbits of particles which cross the horizon- 
tal line between the focal points is cut by the straight 
line from outside, leading indeed again to the selectiv- 
ity in the magic numbers (diamonds in Fig. [5]). Note 
that these calculations were performed with the simple 
reflection program. An analysis similar to that for semi- 
circular hats is also possible for elliptic hats, but shall 
not be carried out here. 



VI. CONCLUSION 

We have been able to show, that mushroom hats with 
openings in the straight boundaries have a very selective 
behavior concerning the number of bounces a particle has 



with the curved boundary before it escapes. Indeed while 
the magic number one always occurs we can to some ex- 
tent design the subsequent numbers and hence the cor- 
responding delay times. We expect that the selectivity 
of the possible delay times in the hat of the mushroom 
should show up not only in classical billiards, but also in 
billiards excited with microwaves. In such experiments 
the open hat mushroom billiard is excited with a single 
pulse input signal and the output time signal is expected 
to be nonperiodic, in contrast to acoustical echoes we 
know from experience and scattering echoes analyzed in 
some detail in Refs. [H, Hl> H3- This opens new av- 
enues in billiard research. An experimental device for 
the measurement of delay times in open mushroom bil- 
liards can be constructed, based on the analogy between 
open mushroom hats and the billiard, with a flat cylin- 
dric microwave resonator of circular shape. The straight 
line defining in the classical billiard the starting and the 
endpoints of orbits which correspond to particles entering 
and leaving the mushroom hat can be realized by a strip 
of microwave absorbing material with an emitting and a 
receiving antenna positioned close to the strip. To test 
how the results presented in this paper carry over to wave 
mechanics, we plan such a microwave experiment with a 
superconducting cavity proceeding along similar lines as 
in the one reported in Ref. (34[ in order to investigate the 
time response. It is though important to repeat, that the 
equivalence between time and bounce numbers holds only 
for small openings. Once we have star like and nonstarlike 
polygon orbits coexisting the corresponding delay times 
get mixed, and while the selectivity of bounce numbers 
may persist, their physical importance may be marginal 
in this case. Also we must be aware, that the analogy 
will break down at long times, but these are probably of 
marginal interest anyway. 
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Appendix: Magic numbers for fixed angular 
momentum 

In this appendix we shall show that for each value of 
angular momentum there exist only three magic numbers 
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n;,n r , and n s where the largest is always the sum of 
the lower ones, n s = ni + n r . For a given length of the 
straight line and a fixed angular momentum, the angles ip 
and x can be calculated from Eq. |T]) and Eq. ^ , where 
< 2ip < ir and < \ < ip [see Eq. ©]. The magic 
numbers are determined by using a procedure based on 
the inequalities Eq. ([7]) and Eq. ((5|) and outlined in the 
main text. 

Now consider Eqs. QJ) and ([6]) which are related to the 
so-called circle map [39| 

<Pn+l = {<t>n + a ) mod 27T . (10) 

For this map it is known for almost 40 years that angles 
4>o originating in a given interval can at most have three 
different recurrence times, i.e., numbers of iterations for 
recurrence to the same interval [39| . The angle parame- 
ter a corresponds to 2ip and <p n to n mod27r, i.e., to the 
orientation of the n-th line segment of the particle orbit, 
where the information on the number of times the parti- 
cle travelled around the circle center until its n-th impact 
with the circle boundary is dropped (see Eq. and the 
remark thereafter). Thus this map indicates one for fixed 
an gula r momentum and we can use the results presented 
in [391 ] to argue that, for a fixed angular momentum, we 
have at most three different delay times. 

An intuitive proof of this result is given here based 
on the translation of the orbit intervals defined in the 
main text after Eq. l[7|). This yields an understanding of 
what happens when the angular momentum varies. As 
an important further result we shall see that for a fixed 
angular momentum less than three delay times occur on 
a subset of measure zero. Since we only consider particle 
orbits which start on the straight line in the O-billiard 
before their first encounter with the circular boundary, 
the first orbit interval coincides with the first line interval 
around zero. The first particles reach the straight line 
when there is the next partial overlap of an orbit interval 
with a line interval at the right hand side or at the left 
hand side. Their particle orbits end there. Accordingly, 
this part of the orbit interval is cut off the orbit interval 
in the next shift. The remaining part evolves until the 
next overlap is reached. 

We first proceed to show in (i) that if at the first over- 
lap after n r reflections of the particle at the circle bound- 
ary, that is shifts of the orbit interval by 2ijj, the orbit 
interval is cut partially from the right hand side and at 
the second overlap after n; reflections from the left hand 
side, then the rest of the interval overlaps with a line in- 
terval at latest after n s = n r + ni reflections. The same 
holds in the opposite case when the orbit interval is first 
cut from the left and then from the right. Next we show 
in (ii) that it is not possible that the orbit interval is 
shortened at the same side for the first two subsequent 
overlaps. As a further step, we demonstrate in (iii) that 
there is no way to find a further overlap between the first 
two overlaps and the one after n s reflections. Finally we 
show, that the possibility of overlapping the complete or- 
bit interval with less than three coincidences of the orbit 



and a line interval is limitted to a subset of measure zero 
within the set of all allowed angular momenta. 

In order to find the magic numbers it is sufficient to 
consider the evolution of the orbit interval center (oc„) 
and the line interval center (lc m ) with increasing numbers 
of shifts n = 0, 1, 2, • • • and m = 0, 1, 2, • • ■ of the orbit and 
line intervals, respectively, where oc n is located at n2ip 
and lc m at mn. The n-th orbit interval partially overlaps 
with the m-th line interval, if the distance between oc n 
and lc m is less than the length of the orbit and the line 
intervals, 2\ [see Eq. (J7J]. It is important to note that 
the length of the orbit and the line intervals is smaller 
than the distance between two subsequent oc's, 2\ < 2ip, 
and the distance is smaller than that of two subsequent 

IC'S, 2lp < 7T. 

(i) Assume that after n r shifts of the orbit interval and 
mi shifts of the line interval the distance between oc Tlr 
and lc mi is smaller than 2%, and that the orbit inter- 
val there partially overlaps at its right end with the line 
interval. Then oc„ r is located to the left of lc mi along 
the angle axis at a distance t\ — min — n r 2ip < 2%. 
The length of the overlapping part of the orbit interval 
equals L\ = 2% — ei- Further assume, that after nj shifts 
the remaining part of the orbit interval partially overlaps 
at its left end with the m,2-th shift of the line interval. 
Then oc„, is located to the right of lc m2 at a distance 
£2 = nflip — rri2ir < 2\ and the length of the overlapping 
part equals L2 = 2\ — (2- The length of the remaining 
part of the orbit interval equals 

L = 2 X -L l -L 2 (11) 
= ei + e 2 -2 X (12) 
= (m — n r )2ip — (mz — mi)ir — 2x- (13) 

If L = 0, then already after two overlaps of an orbit 
with a line interval, all particles will have returned to the 
straight line, such that there are only two magic numbers; 
this case will be treated in part (iv) of the proof. The 
case L < 0, where more than just the remaining part of 
the orbit interval overlaps with the line interval, may be 
excluded, as this would imply that there is another magic 
number n* = (n; — n r ) < ni and an m* = {m% — mi) < 
mi , where the distance between oc„* and lc m * is less than 
2% [see Eq. ([15])]. in contradiction to the assumption. 
For L > 0, the remaining part is adjacent to the line 
interval at its left end; its left endpoint is located at a 
distance x, its right endpoint at a distance t\ + t-2 — x 
from lc m2 . Exactly as after the first n r and mi shifts of 
the orbit and the line interval, after n r and mi further 
shifts of the intervals, the location of oc„ s (n s = n r + 
ni), and therefore that of the remaining part, will have 
moved by a total amount of t\ to the left with respect 
to lc m3 (m.3 = mi + m%) as compared to their initial 
relative locations, i.e., that of oc ni and lc m2 . Accordingly, 
oc„ s is located at a distance |ea — ei| from lc m3 , which 
is smaller than 2\. The left endpoint of the part of the 
orbit interval that remained after n; shifts, is then located 
at a distance \x — £i|, the right endpoint at a distance 
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\ e -2 — x\ from lc m3 . As both distances are smaller than 
X, this part completely overlaps with the line interval. 
Therefore, provided the orbit interval is first cut at its 
right, then at its left end at the magic numbers n r and ni, 
the remaining particles all will have reached the straight 
line again latest after n s — n r + m shifts of the orbit 
interval, or equivalently, bounces of the particle with the 
circle boundary. The same argumentation holds if the 
interval is first cut from the left and then from the right 
hand side. 

If the first magic number equals 1, as in the example 
shown in Fig. dj we already arrived at the end of the 
proof. In this case, the orbit interval will be cut at its 
right end after one shift of the line and the orbit interval 
and oci will be located at a distance ej = tt — 2-0 to the 
left of lc\. At each shift n and m = n of the intervals, the 
distance of oc„ from lc m , which is located to its right, will 
increase by an amount t\ until the orbit interval overlaps 
with the line interval to its left after n; shifts, i.e., until 
the distance between oc„, and lc ni -i is smaller than 2\. 
Hence, at the second overlap the orbit interval will be cut 
at its left end, at a magic number n;. Then, as shown 
in part (i) of the proof, at n s = n; + 1 all particles will 
have reached the straight line, and we obtain the magic 
triplet l,m,ni + 1. The explicit values can be computed 
with Eq. ©. 

(ii) Assume, that the first magic number is larger than 1 
and that the orbit interval overlaps at its left hand side 
after nn and again after n; 2 shifts of oc with respectively 
the rai-th and the m2-th line interval. Then, oc nil and 
oc ni2 will be located to the right of lc mi and lc m2J re- 
spectively, at a distance 

£i = rc a 2V> - rniTT < 2 X (14) 

for the first and 

e 2 = ni 2 2tp - m 2 TT < 2% (15) 

for the second overlap. The length of the first overlapping 
part equals L\ = 2\ — e\, that of the second L 2 = 2x~ e 2 , 
and the orbit interval will be further shortened at the 
second overlap only, if L 2 is larger than L±, that is, if 
e 2 < ei- Then 

< £i - e 2 = (m 2 - rai)7r - (n t2 - nn)2ip < 2\ ■ (16) 

This implies, that there is another magic number n* = 
ni2 — nn < and a number m* = m 2 — mi, where 



the distance ei — e 2 of oc„« from lc m * is less than 2\ 
(right hand side of the inequality) and the orbit interval 
is shortened at its right hand side (left hand side of the 
inequality). This is in contradiction to the assumption. 

(iii) Let the orbit interval be shortened at the right hand 
side after n r reflections and at the left hand side after 
ni reflections. Assume, that the third magic number is 
smaller than the sum of the first and second, n* < n r +ni. 
If after n* shifts the orbit interval is shortened again at 
the left hand side, then, according to (ii), there exists 
an additional magic number n** = n* — n\ < n r , where 
the orbit interval is shortened at the right hand side, in 
contradiction to the assumption. If on the other hand n* 
shortens the orbit interval at the right hand side, then, 
by proceeding as in (ii) , it can be shown that there must 
be another magic number n** = n* — n r < ni where the 
orbit interval overlaps at the left hand side with a line 
interval, again in contradiction to the assumption. 

(iv) Assume that the orbit interval after n r shifts par- 
tially overlaps at its right hand side with a line interval, 
and after ni shifts the remaining part completely overlaps 
with a line interval. This situation corresponds to the 
case L = in part (i), i.e., if (n/ — n r )2ip — (to 2 — toi)7t = 
2\, then there are only two magic numbers. For a given 
length of the straight line, this equality is fulfilled only 
for a set of measure zero from the whole set of allowed 
angular momenta [see Eq. [T]) and Eq. ([4|] . As we vary 
the angular momentum for a fixed length of the straight 
line, the angles ip and x change continuously, and a third 
magic number emerges because the first and second over- 
lap no longer cover the whole orbit interval. There, we 
observe a discontinuity in Fig. [5] 

The arguments given in this appendix only hold for 
finite magic numbers. In the case of the parabolic mani- 
folds described in Ref. [37| only one magic number is left 
finite. The values of angular momentum corresponding 
to those orbits are also of measure zero in the whole set 
of possible angular momenta. In Fig. [SJ we observe a 
singularity at these values of angular momentum. 

Finally we should note, that the proof also holds, if the 
line interval is shifted by an arbitrary angle 2% < f2 < 7r, 
as we only used the fact, that the length of the orbit and 
the line interval is smaller than the distance between two 
subsequent line intervals or orbit intervals. 
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